Free energy and criticality in the nucleon pair breaking process 
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Experimental level densities for "^'^^Yb, ^^'^'^'^^Er, ^'^^'^^'^Dj, and "®'"^Sm are analyzed within 
the microcanonical ensemble. In the even isotopes at excitation energies E < 2 MeV, the Helmholtz 
free energy F signals the transition from zero to two quasiparticles. For E > 2 MeV, the odd and 
even isotopes reveal a surprisingly constant F at a critical temperature Tc ~ 0.5 MeV, indicating 
the continuous melting of nucleon Cooper pairs as function of excitation energy. 
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One of the most spectacular pairing phase transitions 
in nature is the transition from a normal to a supercon- 
ducting phase in large electron systems. The transition 
is triggered at low temperature by massive pairing of two 
and two electrons into spin ,7 = pairs, so called Cooper 
pairs [1]. 

For atomic nuclei, the pairing phase transition is ex- 
pected to behave differently. First of all, the nucleus is an 
isolated, few body system with two species of fermions. 
Surface effects are prominent, and the coherence length 
of nucleons coupled in Cooper pairs is larger than the 
nuclear diameter. Furthermore, there are non-negligible 
energy spacings between the single particle orbitals. All 
these facts make the nucleus an inherently small system. 
Also, other types of residual interactions than pairing are 
of importance. The influence of these peculiar constraints 
on the nucleus has been investigated theoretically for a 
long time [2-5], however, only limited experimental in- 
formation is available to describe the nature of pairing 
within the nucleus. 

The Oslo group has developed a method to derive si- 
multaneously the level density and 7-ray strength func- 
tion from a set of primary 7-ray spectra [6] . The method 
has been well tested and today a consistent data set for 
eight rare earth nuclei is available. In the present work 
we report for the first time on a comprehensive analy- 
sis of the evolution of the pairing phase transition as a 
function of the nuclear excitation energy. 



II. EXPERIMENTAL LEVEL DENSITIES 



Level densities for ^^^^'^'^^Yh, 166,167^^.^ I6i,i62j3y^ 
I48,i49gj-|^ have been extracted from particle-7 coinci- 
dences. The experiments were carried out with 45 MeV 
■^He projectiles accelerated by the MC-35 cyclotron at 
the University of Oslo. The data were recorded with the 
CACTUS multidetector array using the pick up ('^IIe,Q;7) 
reaction on i^^.tTSyb, i67Er, I62,i63j3y^ i49g^ 

gets and the inelastic (^He,^He'7) reaction on ^^^Er and 
^■*^Sm targets. The charged ejectiles were detected with 
eight AE-E particle telescopes placed at an angle of 45° 
relative to the beam direction. Each telescope comprises 
one Si front and one Si (Li) back detector with thicknesses 
of 140 and 3000 /im, respectively. An array of 28 Nal 7- 
ray detectors with a total efficiency of ~15% surrounds 
the target and particle detectors. From the reaction kine- 
matics, the measured ejectile energy can be transformed 
into excitation energy E. Thus, each coincident 7 ray 
can be assigned to a 7 cascade originating from a specific 
energy E. These spectra are the basis for the extraction 
of level density and 7-strength function as described in 
Ref. [6]. Several interesting applications of the method 
have been demonstrated, see, e.g., Refs. [7-10]. 

The level densities for i^^'^^yb, 166,167^^.^ I6i,i62j3y^ 
and ^^^'^"'^Sm are shown in Fig. 1. The level densities are 
normalized at low excitation energies where (almost) all 
levels are known, and at the neutron binding energy _B„ 
where the level density can be estimated from neutron- 
resonance spacings. The spin window populated in the 
reactions is typically / ~ 2-6h. Already, three general 
comments can be made to these data: («) above 2 MeV 
excitation energy, all level densities are very linear in 
a log plot, suggesting a so-called constant temperature 
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level density, (ii) the level densities of the odd-even iso- 
topes are larger than for their neighboring even-even iso- 
topes, and (iii) the even-even isotopes show a strong in- 
crease in level density between 1 and 2 MeV, indicating 
the breaking of Cooper pairs. 

It should be noted that the transitions considered 
here are low temperature phenomena. The ^^^-^^^Yb, 
I66,i67gj.^ and i6i,i62Qy nuclei have well deformed shapes, 
and various calculations in this mass region [11 14] indi- 
cate that the transition from deformed to spherical shape 
occurs at much higher temperatures than the tempera- 
ture at which the first pairs break. However, for nuclei 
closer to the = 82 shell gap, e.g., ^^^'^'''^Sm, the co- 
existence between deformed and spherical shapes at low 
temperatures cannot be excluded, as discussed in Ref. 
[15]. 



III. FREE ENERGY AND CRITICAL 
TEMPERATURE 

The statistical microcanonical ensemble is an appro- 
priate working frame for describing an isolated system 
like the nucleus. In this ensemble the excitation energy 
E is fixed, in accordance with the observables of our ex- 
periments. The microcanonical entropy is given by the 
number of levels f2 at i? 

S{E) = kBlnn{E), (1) 

where the multiplicity Q is directly proportional to the 
level density p by Q{E) = p[E)/ pQ. The normalization 
denominator pQ is adjusted to give 5 ~ for T ~ 0, 
which fulfills the third law of thermodynamics. Here, we 
assume that the lowest levels of the ground state bands of 
the ^'^^Yb, ^'^^Er, ^^^Dy, and "^Sm nuclei have tempera- 
tures close to zero, giving on the average po = 2.2 McV^^ . 
In the following, this value is used for all eight nuclei and 
Boltzmann's constant is set to unity [ks = !)• 

In order to analyze the criticality of low temperature 
transitions, we investigate the probability -P of a system 
at the fixed temperature T to have the excitation energy 
E, i.e., 

P{E, T, L) = n{E) exp {-E/T) /Z{T), (2) 

where the canonical partition function is given by Z{T) — 
n{E')exp{-E'/T) dE'. Implicitly, the multiplicity 
of states ^l{E) depends on the size of the system, denoted 
by L. Often, it is more practical to use the negative 
logarithm of this probability A{E,T) = - In P{E,T), 
where we in the following omit the L parameter. Lee 
and Kosterlitz showed [16,17] that the fimction A{E,T), 
for a fixed temperature T in the vicinity of a critical 
temperature Tc of a structural transition, will exhibit 
a characteristic double-minimum structure at energies 
El and i?2- For the critical temperature T^, one finds 
A{Ei,Tc) = A{E2,Tc). It can be easily shown that A is 



closely connected to the Helmholtz free energy and the 
previous condition is equivalent to 

F,{Ei) = F,{E2), (3) 

a condition which can be evaluated directly from our ex- 
perimental data. Here, it should be emphasized that Fc 

is a linearized approximation to the Helmholtz free en- 
ergy at the critical temperature Tc according to Fc{E) = 
E—TcS{E), thereby avoiding the introduction of a caloric 
curve T{E). The free-energy barrier at the intermediate 
energy Em between Ei and E2, is given by 

AF, = Fc{E„,) - Fc{Ei). (4) 

Now, the evolution of Ai^c with increasing system size L 
may determine the order of a possible phase transition 
[16,17]. These ideas have, e.g., recently been applied to 
analyze phase transitions in a schematic pairing model 

Figure 2 displays a schematic description of the en- 
tropy for even-even, odd-mass and odd-odd nuclei as 
function of excitation energy. In the lower excitation en- 
ergy region of the even-even nucleus, only the ground 
state is present, and above E ^ 2A the level density is 
assumed to follow a constant temperature formula. It 
has been shown [19,20] that the single particle entropy s 
is an approximately extensive thermodynamical quantity 
in nuclei at these temperatures. The increase in entropy 
at the breaking of the first proton or neutron pair, i.e., at 
E = 2A, is roughly 2s in total for the two newly created 
unpaired nucleons. The requirement Fc{Ei) = Fc{E2), 
where T = T^, gives E2-E1 = [S{E2) - S{Ei)]. Thus, 
with the assumed estimates above, we obtain the rela- 
tion A = Tc-s, which may be used to extract the critical 
temperature for the pairing transition. Adopting typical 
values of A = 1 McV and s — 2fcs [19,20], we obtain 
Tc = 0.5 MeV. For the odd-mass case, one starts out 
with one quasiparticle which gives roughly one unit of 
single-particle entropy s around the gromid state. The 
three quasiparticle regime appears roughly at = 2A 
with a total entropy of 3s. The region between E = 
and 2 McV is modeled with a step at ^ 2 McV, how- 
ever, in real nuclei the level density is almost linear in a 
log-plot for the whole excitation energy region due to the 
smearing effects of the valence nuclcon. In the case of an 
odd-odd nucleus, one starts out with two units of single- 
particle entropy. The two valence nucleons represent a 
strong smearing effect on the level density and the mod- 
eled step striicturc in entropy at E = 2A for the onset of 
the four quasiparticle regime is completely washed out. 

In the higher excitation region, further steps for transi- 
tions to higher quasiparticle regimes are also washed out 
due to the strong smearing effects of the already present, 
unpaired nucleons. The slope of the entropy with excita- 
tion energy is determined by two competing effects: the 
quenching of pairing correlations which drives the cost 
in energy lower for the breaking of additional pairs, and 
the Pauli blocking which reduces the entropy created per 
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additional broken pair. The competing influence of both 
effects is modeled by a constant-temperature level den- 
sity with the same slope for all three nuclear systems and 
a slightly higher critical temperature. In this region of 
the model, there are infinitely many excitation energies 
where the relation Fc{Ei) = Fc(i?2) is fulfilled. 

The breaking of proton or neutron pairs are thought 
to take place at similar excitation energies due to the 
approximate isospin symmetry of the strong interaction. 
It is indeed commonly believed that the pairing gap pa- 
rameter A and, thus, the critical temperature Tc for the 
breaking of Cooper pairs, are approximately the same 
for protons and neutrons. Furthermore, interactions be- 
tween protons and neutrons will certainly wash out any 
differences in behavior between the proton and neutron 
fluids. In Fig. 3, the influence of differences in pro- 
ton and neutron pair breaking is investigated within our 
schematic model. Here, we assume that neutrons break 
up at 2A creating an entropy of 2s. The protons are as- 
sumed to break up at 10% higher excitation energy (since 
Z < N) creating 10% less entropy (due to the larger pro- 
ton single particle level spacing). The entropy of the total 
system of either proton or neutron pair breaking gives 



S{E)=\n [e^''(^)+e^''(^)l -ln2, 



(5) 



where the last term assures that S" = in the ground 
state band. The requirement Fc{Ei) = Fc{E2) gives 
ri"^ = A/s for neutrons (as in Fig. 2) and t]^^ = 
l.lA/0.9s = 1.22A/s for protons. In the combined 
system of both neutrons and protons, a value of = 

(Ti"^ -hTi^V2 = l.llA/s is deduced from Fig. 3. Thus, 
typical fluctuations in the pairing gap parameter and the 
single particle entropy for neutrons and protons give only 
small changes in the extracted critical temperature. 



IV. EXPERIMENTAL RESULTS 

In order to experimentally investigate the behavior for 
even isotopes, linearized free energies Fc for certain tem- 
peratures Tc are displayed in Fig. 4. The data clearly 
reveal two minima with Fc{Ei) = Fc{E2) = Fq, which 
is due to the general increase in level density around 
E 1^ 2 MeV, as schematically shown in Fig. 2. For all 
nuclei, we obtain Ei ^ and E2 ^ 2 MeV which com- 
pares well with 2A. We interpret the results of Fig. 4 as 
the transition due to the breaking of the very first nucleon 
pairs. The deduced critical temperatures are Tc — 0.47, 
0.40, 0.47 and 0.45 MeV for ^^^yb, 166^^.^ i62jQy^ 
"^^^Sm, respectively. 

Recently [7] , another method was introduced to deter- 
mine the critical temperatures in the canonical ensemble. 
Here, the constant temperature level density formula for 
the canonical heat capacity Cv(T) = (1 — T/r)"^ was fit- 
ted to the data in the temperature region of 0-0.4 MeV 
corresponding to excitation energies between 0-2 MeV, 



and the fitted temperature parameter t was then iden- 
tified with the critical temperature Tc- Since a constant 
temperature level density formula implies a constant lin- 
earized Helmholtz free energy Fc{E) (provided r = Tc), 
this former method is almost equivalent to the present 
method, i.e., of identifying the temperature Tc for which 
the linearized Helmholtz free energy is in average con- 
stant. Therefore, it is not surprising that the extracted 
critical temperatures Tc = 0.49, 0.44, 0.49 and 0.45 MeV 
for the respective nuclei using the older method [7,8,10] 
coincide well with the critical temperatures presented 
in this work. However, while the previous method was 
based on an ad-hoc assumption of the applicability of a 
constant temperature level density formula, the present 
method has a much firmer theoretical foundation. 

The height of the free-energy barrier should show a dif- 
ferent dependence on the system size L according to the 
order of a possible phase transition [16,17]. The barriers 
deduced from Fig. 4 yield AFc ~ 0.5-0.6 MeV, values 
which seem not to have any systematic dependence on 
the mass number A within the experimental uncertain- 
ties. Even with better data, an unambiguous dependence 
of the barrier height on the system size would be unlikely 
when using ^ as a measure for the parameter L since the 
relevant system size for the very first breaking of Cooper 
pairs might be characterized by only a few valence nu- 
cleons. Another complicating interference is that other 
properties of the nuclear systc;ni which might influence 
the onset of pair breaking also change with mass num- 
ber, e.g., deformation, pairing gap, and locations of single 
particle levels around the Fermi surface. 

In the schematic model of Fig. 2, we would expect 
a free energy barrier of AF^. = 2A <^ 2 MeV E = 
2A ^ 2 MeV. However, the data are more consistent 
with the dotted lines of Fig. 2 indicating a smoother be- 
havior around the expected steps due to the existence of 
collective excitations like rotation and /3, 7, and octupole 
vibrations between 1 and 2 MeV for the even-even nuclei, 
and due to the increasing availability of single-particle 
orbitals for the odd nucleon in the case of odd nuclei. 
Thus, we expect the centroid of the barrier to be shifted 
down in energy with a corresponding proportional re- 
duction of the barrier height, and an inspection of Fig. 4 
indeed shows that the free energy barrier is 0.5-0.6 MeV 
at ~ 0.6 MeV excitation energy for the even-even nuclei. 
A similar analysis of the odd isotopes is difficult to ac- 
complish since there seems not to be any common struc- 
tures. Here, the unpaired valence neutron smears out the 
effects of the depairing process too much to be visible in 
the present data. However, it has been attempted in Ref. 
[8] to interpret the structure in the level density of ^^''Er 
around 1 MeV in terms of a first order phase transition. 

The smearing effect is expected to be even more pro- 
nounced for the breaking of additional pairs. Figure 5 
shows the linearized Helmholtz free energy for all eight 
nuclei investigated, but at slightly higher critical tem- 
peratures than in Fig. 4. The critical temperature Tc 
is found by a least fit of a constant value Fq to the 
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experimental data. The fit region is from E = 2 MeV 
up to 5 and 7 MeV for the odd and even isotopes, re- 
spc;ctively, giving normalized values in the range from 
0.5 to 2.5. Here, instead of a double-minimum struc- 
ture, a continuous "minimum" of Fc is displayed for sev- 
eral MeV. This observation allows us to conclude that 
the further depairing process cannot under any circum- 
stances be interpreted as an abrupt structural change in 
the nucleus typical for a first order phase transition. The 
constant lines of Fig. 5 visualize how surprisingly well 
Fq fits the data: the deviations are typically less than 
100 keV.""^ The ongoing breaking of further Cooper pairs 
overlapping in excitation energies above 2 MeV is there- 
fore contrary to what is found in the schematic model of 
Ref. [18] . This is probably due to strong residual interac- 
tions in real nuclei, like the quadrupole-quadrupole inter- 
action, which were not taken into account in the model 
calculation. Thus, (nearly) all excitation energies above 
2 MeV will energetically match with the costs of break- 
ing nucleon pairs. Here, all excess energy goes to the 
process of breaking pairs. Since the gain in entropy dS 
is proportional to dE, the microcanonical temperature, 
T{E) = {dS/dE)~^, remains constant as function of ex- 
citation energy, and the level density displays a straight 
section in the log plot. 

At higher excitation energies than measured here, the 
pairing correlations vanish and the system behaves more 
like a Fermi gas. Here, the free energy will indicate the 
closing stage of the depairing process by increasing F,,. 
with Fc> Fq. However, in this regime also shape transi- 
tions and fluctuations as well as the melting of the shell 
structure may play a role and give deviations from a 
simple Fermi gas model with p cx exp(2-\/ aE), a being 
the level density parameter. Unfortunately, these very 
interesting phenomena cannot be investigated with the 
present experimental data. 

The fitted value Fq contains information on the en- 
tropy of the system at Tc through S = [E — Fq)/Tc. 
In Fig. 6, we have compared the entropy for the various 
nuclei at an excitation energy E = 4: MeV, an energy 
where all nuclei seem to "behave" equally well (see Fig. 
5). Figure 6 also shows that the odd mass nuclei display 
generally higher entropy regardless of the mass number 
A being one higher or lower than the neighboring even 
isotope. We also observe that since the ^'*^'^^^Sm nuclei 
are not niid-shcll nuclei, they show less entropy, reflect- 
ing the lower single particle level density when approach- 
ing the TV = 82 shell gap. By evaluating the odd-even 
5S = Sodd — S'ovcn, we find SS ^ 2 for all four isotopes, as 
shown in the lower panel of Fig. 6. This means that ex- 
cited holes and particles have the same degree of freedom 
with respect to the even mass nuclei. 



V. CONCLUSION 

Unique experimental information on level densities for 
eight rare earth nuclei is utilized to extract thermody- 
namic quantities in the microcanonical ensemble. The 
linearized Helmholtz free energy is used to obtain the 
critical temperatures of the depairing process. For a crit- 
ical temperature just below Tc ^ 0.5 MeV, we observe a 
structural transition of even nuclei in the E = — 2 MeV 
region due to the breaking of the first nucleon pair. Un- 
fortunately, it was not possible to use the development 
of the barrier height AFc with the size of the system L 
to conclude on the presence of a thermodynamical phase 
transition and its order. The critical temperature for the 
melting of other pairs is found at slightly higher temper- 
atures. Here, we obtain a surprisingly constant value for 
the linearized Helmholtz free energy, indicating a con- 
tinuous melting of nucleon Cooper pairs as function of 
excitation energy. The conspicuous absence of a double- 
minimum structure in Fc for this process is at variance 
with the presence of a first-order phase transition in the 
thermodynamical sense. The entropy difference between 
odd and even systems is found to be constant with respect 
to excitation energy and is consistent with the expected 
values of the single particle entropy in these nuclei. 
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FIG. 1. Experimental level densities for the nuclei ^^^'^^Yb, ^^^'^^^Er, ^^^^^^^By, and "^■"^Sm. The data are taken from 
Refs. [7,8,10]. 
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FIG. 2. Schematic representation of the entropy S in units of the single particle entropy s (top panel) for even-even (solid 

line) odd-mass (dash-dotted line) and odd-odd (dashed line) nuclei. For the purpose of the figure, the steps in entropy are 
drawn slightly staggered in energy. Lower panel: linearized Hclmlioltz free energy Fc at the critical temperature Tc of even-even, 
odd-mass and odd-odd nuclei. All energies are measured in units of the pairing gap parameter A. The dotted lines indicate 
the situation if additional levels are included below the steps in entropy. 
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FIG. 3. Same as previous figure in the case of an even-even nucleus but for unequal proton and neutron fluids. Curves are 
given for the neutron fluid alone (dashed lines with pairing gap parameter A and single particle entropy s), the proton fluid 
alone (dotted lines with 1.1 A and 0.9s) and the composite system (solid lines). 
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FIG. 4. Linearized Helmholtz free energy at the critical temperature Tc. The constant level Fo connecting the two minima 
is indicated by lines. 
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FIG. 5. Linearized Helmholtz free energy at the critical temperature Tc. The fitted constant level Fo is indicated by lines. 
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FIG. 6. Experimental entropy evaluated in the microcanonical ensemble at excitation energy E = A MeV and temperature 
Tc. In the lower panel the odd-even difference 5S = displayed for the four isotopes. 
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